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Abstract 

TJ . ... 

We examine first-order perturbative results based on jellium for the surface energy of 
slabs of simple metals, using various local pseudopotentials (Ashcroft, Heine-Abarenkov 
and evanescent core). The difference between the pseudopotential and the jellium potential 
is averaged along the plane parallel to the surface. We compare these perturbative results 
with those of the stabilized jellium model (a modification of the regular jellium model in 
which the perturbation appears in the energy functional right from the outset) and with 
the output of other perturbative and non-perturbative calculations. 
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1 Introduction 

The jellium model is the simplest model which can be used to describe simple metals (metals 
with s or p valence electrons). It avoids the atomic non- uniformities by replacing the ionic cores 
by a positive uniform background. It describes qualitatively the work function, but it predicts 
negative surface energies for metals with high valence-electron density. 

Lang and Kohn, who were the first to apply the jellium model to surfaces 012], introduced 
a perturbative correction to get realistic results. This correction was simple enough to keep most 
of the original simplicity: the self-consistent density was still that of jellium and the perturbation 
(difference between the lattice potential and that of the uniforme positive background) was av- 
eraged over the surface plane. They also improved on the description of the ion-ion interaction 
through the so-called classic cleavage surface energy. In this way, they obtained face- dependent 
surface energies, which were always positive. 



Later, the influence of the discrete ions in the electronic density was taken into account [3] 
via second-order perturbation theory, which includes the linear response of the electronic distri- 
bution to the lattice potential [I]. Rose and Dobson [2] were the first to work out second-order 
surface energy terms, but they used the linear response of bulk jellium, in a kind of local density 
approximation. Second-order perturbation theory using the linear response of a jellium slab in 
the Random Phase Approximation has been worked out by Barnett and coworkers [H] and also 
by Eguiluz jE]. These calculations, which are three-dimensional, show a noticeable influence of 
the second-order term in the face-dependent surface energies. While the surface energies depend 
strongly on the exposed face in the first-order perturbative model of Lang and Kohn, the second- 
order results show a weaker dependence. 

The stabilized-jellium model or structureless pseudopotential model [3IHIE]; a modification 
of the regular jellium model in which the perturbation appears in the energy functional right from 
the outset, includes the perturbation in the effective potential of the self-consistent Kohn-Sham 
equations and, as a consequence, its effect in the electronic density. Originally, this model was 
intended to describe flat surfaces and therefore did not include any structure in the averaged 
perturbation. However, the dependence of the surface energy on the atomic corrugation of a 
particular face of a metal was incorporated multiplying the flat surface results by a term based 
on the liquid-drop model [3 QT1] . The results showed a much weaker face dependence than in the 
work of Lang and Kohn, in agreement with second-order perturbative results. Considering a dipole 
barrier, due to corrugation, the model was also adapted to improve the previous face- dependent 
surface energies [TT] and to calculate face- dependent work functions GDI HI] • 

More elaborated, and therefore more computationally demanding non-perturbative calcu- 
lations, are now available for surfaces j!2j . The experimental difficulties to get surface energies of 
the different surfaces to compare with the perturbative results based on jellium were then obviated 
by the predictions of these full atomistic calculations. The results of the stabilized-jellium model 
were found to be fairly realistic for several metals [13J. 

Some of the above-mentioned calculations were performed for slabs or thin films, i.e., sys- 
tems made out of a few atomic layers. Slabs, which are convenient to obtain second-order per- 
turbative results for surfaces, are interesting in their own since they exhibit quantum size and 
self-compression effects. Jellium slabs, showing quantum size effects, were examined in a seminal 
paper by Schulte [Hj. The second-order perturbative energies of Barnett et al. |Hj and Eguiluz [Oj 
were obtained for slabs. Slabs of stabilized jellium were recently examined by us [13]. We have 
shown that they are able to describe both quantum-size and self- compression effects. 

The above perturbative treatments were implemented using Ashcroft's empty-core pseu- 
dopotential [IE] or the local Heine- Abarenkov pseudopotential [Uj. In this paper we study alu- 
minum slabs, with 9 and 17 layers, along the lines of Lang-Kohn's perturbation theory using the 
Ashcroft, the Heine- Abarenkov, and the recent evanescent- core pseudopotential |18j . which has 
the advantage of having a smooth repulsion. We compare our results with those of stabilized 
jellium and with other perturbative and non-perturbative calculations. Finally, we refer to the 
possibility of considering stabilized jellium as a zero-order description of a metal surface. 



The ultimate goal of our research is to perform systematic first and second-order pertur- 
bative calculations for surfaces of simple metals with the evanescent-core pseudopotential. This 
will be an extension to surfaces of the systematic treatment of the energetics and mechanical 
properties we have made for metallic solids in different crystal structures [19J. Although ab initio 
calculations are nowadays clearly the method of choice for bulk or surface systems, perturbative 
treatments still have their role to describe trends along the the periodic table and along different 
crystallographic structures. Above all, they can provide understanding of the physics of metal 
cohesion. 

The outline of this paper is as follows. In Section 2, the theoretical background is provided, 
in Section 3 the perturbative results are presented. Conclusions appear in Section 4. 



2 Perturbative corrections to jellium model 

2.1 Electronic and Madelung subsystems 

In order to study a neutral metallic system formed by a fixed lattice of ions and electrons 
with density n{f), we start with a superposition of two simple systems: the valence electrons 
moving in a positive background of density n + (r) = n inside the metal (n + (f) = outside), where 
n is the mean electronic density, and the ions lattice embedded in a negative background with the 
same density and size as the positive one. 

The first of these subsystems, referred to as the electronic subsystem, is the jellium model. 
We denote its ground-state energy by Ej[tiq\, n being its electronic density which we assume to 
be a reasonable approximation to n(f). 

The second subsystem, called Madelung subsystem, has energy Em- If one wishes to correct 
the jellium description perturbatively, one should consider the interactions within the Madelung 
system and between this and the electronic subsystem, Em,b- The total energy is 

E = E J [n Q ]+E M + E Mfi . (2.1) 

The energy of the jellium model (electronic subsystem) is, in atomic units, given by the 
following density functional: 
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where T[no] is the noninteracting kinetic energy, i^cbio] the exchange-correlation energy (often 
evaluated in the local density approximation, which is incidentally adequate for jellium surfaces 



J20JI2I]), an d the other terms represent, respectively, the electronic repulsion, the self-repulsion of 
the positive background and the attraction between the electrons and the positive background. 
On the other hand, the Madelung energy is given by 
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where the summations run over the ionic positions R(l) and Z is the charge of each ion. The first 
term is the self-repulsion of the negative background, the second the Coulomb repulsion of the 
ions and the third the interaction between the negative background and the ions. 



2.2 First-order perturbative correction 

In order to simplify the calculations, it is convenient to use a local pseudopotential, v ps , 
to represent the ion-electron interactions (from now on, by electrons we mean valence-electrons). 
The pseudopotential, due to a pseudo-ion of charge Z located at the position R(l), can be written 
as a sum of two contributions: an attractive long-range Coulomb part and a short-range repulsive 
part: 

Z 
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The energy arising from the interactions between the Madelung and the electronic subsys- 
tems can be written in first order as: 
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The terms in the first brackets of (|2.5j) represent the interaction between the Madelung 
subsystem and the electrons of the electronic subsystem. They correct the jellium-electron inter- 
actions (last term in the right-hand side of (J2.2J) ). The difference 5v (r) between the potential of 
a pseudo-ions lattice and the potential of the jellium background appears in the first-order term, 
which is the first correction to the jellium model: 



ES?[n Q ]= d 3 r5v(r)n (r), 
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The second brackets of ()2.5|) include terms representing the interaction between the 
Madelung subsystem and the positive background of the electronic subsystem. Taking advan- 
tage of the pseudopotential form (J2.4J) we may rewrite it as 



E M)+J = - Jd 3 r5v(r)n + (r) + f d 3 r^R (f- R(lj)n+(r). (2.8) 

If we add E M of (|2.3JI to E M+J (as written in the second brackets of 1)2.50 ). we find the 
second correction to the jellium model, the one which improves the description of the ion-ion 
interaction: 
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which amounts to replacing the self-repulsion of the background by the repulsion of the point ions. 
The described perturbative approach is based on general density functionals and may be 
applied as well to solids, surfaces, slabs or clusters. Lang and Kohn PQ |2] used it to calculate 
surface properties of metals. However, they replaced the difference potential Sv(r) in ()2.6|) by 
its average along the plane parallel to the surface, Sv(z) (z being the direction perpendicular to 
the surface), keeping therefore the one-dimensionality of the underlying jellium model, for which 
rio(r) = n (z). They wrote: 

8v(z) = (5v(r)) = J2 n 2nd I ' dx\\ v ps (\x\ + \z - ^(A;)| 2 ]" 1/2 )x|| - <p + (z), (2.10) 

k J 

where k runs over the ionic planes, d denotes the inter-layer spacing, xu is the distance in the 
plane parallel to the surface from any ion, R z (k) is the position of the kth plane of ions, and (p+(z) 
is the potential due to the uniform positive background. 

Lang and Kohn used the Ashcroft empty-core pseudopotential ^E] > which requires for each 
metal a single parameter, the core radius re'- 
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A better pseudopotential is the local form of the Heine- Abarenkov pseudopotential [T7j . 
which includes two parameters, Re and u, the first measuring the core radius and the second the 
amount of repulsion in the core, 

ps V U V \- Z /\f-R(l)l \f-R(l)\>R c . y } 

A recent pseudopotential incorporating an exponential decay of the core repulsion, devised 
for evaluating systematically the energetics of simple metals, is the evanescent core pseudopoten- 



tial [TBI, which also depends on two parameters, R and a: 
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where A and .B are simple functions of a. The smoothness of this potential assures good con- 
vergence of its Fourier transform and its suitability to second-order perturbative calculations. It 
yields overall good results for simple metal solids and clusters ^Bl 033 1221 I2HJ , when the parameters 
R and a are fitted to solid-state information. 

Perturbative first-order theory can be applied to slabs along the lines of Lang and Kohn 
using these or other local pseudopotentials. It can also be done exactly, i.e., without taking any 
average [SI E], but then part of the simplicity of the jellium model is lost. 

If one is working with slabs, the surface energy can be extracted from the total slab energy 
per unit area by subtracting the corresponding bulk energy: 

a(L) = ^ A [E(L)-nLAe bulk ], (2.14) 

where A is an area, L is the slab width (the width of the jellium background in our system) and 
e buik j g ^Q total energy per particle in the bulk. We can decompose the slab surface energy in 
various parts. For the functional of (|2.1J) and considering (J2.5J) . the surface energy reads as: 



a = aj + cr M + a^J + a M ,+j- (2.15) 

The first term in the right-hand side is the jellium surface energy. The second is the Madelung 
surface energy, for which the following classical cleavage formulae may be used, as Lang and Kohn 
did 

<jm ~ aZn, (2-16) 

where a is a tabulated constant for each face of a given crystal structure. The third term is the 

first-order pseudopotential surface energy, which may be approximated using the average value 

for the perturbative potential given by (J2.10J1 . 

1 f°° 
cr« ~ -J dz 6v(z)[n(z) - n + (f)]. (2.17) 

Finally, <Jm,+j is a cleavage piece which is different from zero only when the ion cores appear 
out of the jellium surface. We may use, for this term, the expression of [5], where the potential is 
averaged in planes parallel to the surface, as previously done for the difference potential: 

a M}+J ~-±J dzfeu R (r-R(lj)\n, (2.18) 

where the angular brackets denote the surface average. 

When the width of the slab approaches infinity, the surface energy of the slab approaches 
the surface energy of the semi-infinite system. 



2.3 Stabilized jellium 

In the stabilized-jellium model, which is based on the perturbative-variational concept of 
Monnier and Perdew j^j, the perturbative potential, conveniently averaged, is included in the 
Kohn-Sham equations. The new self-consistent density is a better approach to the real density 
than the jellium one. The corrections to the jellium description of the ion-ion and electron-ion 
interactions are now averaged out: the self-repulsion of the jellium positive background inside 
each Wigner-Seitz sphere e is subtracted (the jellium is supposed to describe well the repulsion 
between cells) and, in the perturbative energy, Sv(r) is taken to be constant inside the metal and 
equal to its average over the volume of the Wigner-Seitz sphere, {Sv) ws . The stabilized jellium 
energy functional is given by 

E S j[n] =Ej[n]-e d 3 r n+(f) + (Sv) ws / d 3 r 6(r)n(r) 

or, using the equality [7j e = (Sv)ws -%- u>r, where cm is the (bulk) Madelung energy per 
particle and wr is the pseudopotential repulsion averaged in the Wigner-Seitz sphere, 

E SJ [n] = Ej[n] + (e M + w R ) f d 3 rn + (f) + (5v) ws f d 3 r ^±p [n(r) - n+(r)] . (2.19) 

While in the jellium model the energy per particle of the bulk system has a single minimum 
at a density close to that of sodium 

J bulk 

dn 

in the stabilized-jellium model that energy has a minimum for each metal at the corresponding 
experimental density 

J c bulk 

de sJ I =0 

j \n=n exp u. 

This condition is fulfilled by adjusting a pseudopotential parameter. 

The second term in the right-hand side of (|2.19|) does not contribute to the surface energy. 
The stabilized jellium functional has been recently applied to metal slabs [IHj, the conclusion being 
that the surface energy obtained from (J2.14J) and (J2.19J) 

o~sj = cr S j,j + cr S j,ps (2.20) 

gives a reasonable description of aluminum slabs in comparison with ab initio results, but it fails 
for lithium slabs, a case where the necessity for a non-local pseudopotential is known and for 
which first-principles calculations showed untypical features for a simple metal |24j . Note that 
the first term of (20) is similar to the jellium surface energy but is evaluated with the self-consistent 
density obtained from the functional (|2.19J) . On the other hand, the second term differs from the 
Lang-Kohn perturbative term by the use of a 3D average for the perturbative potential instead 
of a 2D one and by the use of the self-consistent stabilized jellium electron density instead of the 



jellium electron density. In summary: although inspired by it, the stabilized jellium model is not 
perturbation theory. 

In ^H], we have shown that the application of the stabilized jellium model to slabs leads to 
quantum size effects, i.e., fluctuations in the surface energy and work function, which are similar 
to those known for the jellium model but are around more realistic values. Moreover, we have 
shown that, fixing the width of the slab, energy minimization with respect to background variation 
leads to a higher background density inside the slab, i.e., the system tends to self-compress. 

Using functional (|2.19|) no difference shows up between different crystallographic faces. 
However, ad-hoc modifications of stabilized jellium have been proposed to describe the difference 
between various exposed faces EH E]. In these approaches, the self-consistent density is 
obtained by considering a face-dependent potential but, for the sake of realism, the use of the 
latter is avoided in the final evaluation of the first-order surface energy. This methodology describes 
reasonably well the face-dependence found by more sophisticated methods. 

A modification of stabilized jellium has been made by Montag, Reinhard, and Meyer [2*K] . 
who tried to incorporate the cleavage energy in a phenomenological way (by fitting to empirical 
surface energies). 



3 Results 

We have considered jellium slabs corresponding to 7 and 19 layers of aluminum (fee lattice), 
cut along the three main planes (111), (100) and (110) (by decreasing order of planar density and, 
therefore, by increasing order of interplanar distance). We expect perturbation theory based on 
jellium to converge better for the planes which are most close-packed and therefore more similar 
to a flat surface. 

Fig. 1 represents our three pseudopotentials for aluminum. The following values for the 
pseudopotentials parameters have been used (all in atomic units, except a, which is dimension- 
less): r c =1.12 (Ashcroft); R C =1A017 and u=-0.3921 (Heine- Abarenkov); and #=0.317, a=3.512 
(evanescent core). The first value is simply derived from a stability condition for the bulk energy 
within first-order perturbation theory. The second pair of values come from a stability condi- 
tion for the bulk energy at the experimental density and from matching the bulk modulus (within 
second-order perturbation theory) to the experimental value jH] • The last pair of values arise again 
from a stability condition for the bulk energy within second-order and from matching first-order 
values [THJ to all-electron values for the number of electrons in the interstitial region (zone between 
the Wigner-Seitz cell and the inscribed sphere). It has been shown that the latter requirement 
does not differ much from the demand for a realistic bulk modulus JIH] • 

Jellium surface energies Gj are straightforward to evaluate. The cleavage energy gm has 
been taken from jU] . We note that the differences between the values of gm for various faces are 
large. The term gm,+j is small: it may be even zero depending on the size of the core radius [9 j . 
The perturbative potentials Sv(z) corresponding to the different pseudopotentials are represented 



in Fig. 2. They enter in the calculation of a\, s . We note the better smoothness of the Sv(z) arising 
from the evanescent core potential. 

Table 1 shows total surface energies, together with their components, obtained with the 
three pseudopotentials for the slabs with 7 and 19 layers. The table illustrates the importance of 
the positive Madelung contribution. For the same pseudopotential the surface energy increases 
when going from the (111) face to the (110) face. The results show a strong dependence on 
the pseudopotential, with the Ashcroft result, which in principle is the most unrealistic, being 
discrepant from the other two (it is always bigger). Without the Madelung term, the perturbative 
correction of Lang and Kohn would lead to a positive surface energy only for the (111) face, 
precisely that considered by those authors. 

Table 2 shows the flat stabilized jellium results for aluminum slabs with thicknesses 
corresponding to 7 and 19 layers. The total surface energies are in general very different not only 
from the jellium ones but also from first-order perturbation results. 

In spite of the quantum size oscillations, slabs may be used to estimate the surface en- 
ergy of the semi-infinite system. Indeed, the slab with thickness corresponding to 19 layers is a 
good approximation to the latter. Table 3 allows for comparing our previous slab results with 
semi-infinite results (Lang-Kohn, stabilized jellium, and other perturbative and non-perturbative 
results). Actually, the non-perturbative result was obtained with the plane- wave pseudopoten- 
tial method for 12 layers separated by a vacuum of 6 layers, but this should represent well the 
semi- infinite system. Accepting this non-perturbative as the best result, we have to conclude that 
its agreement with the (111) surface energy of Lang and Kohn, who used the Aschroft potential, 
is accidental. The same pseudopotential for any other face leads to disagreement. However, the 
most striking conclusion of Table 3 is that the face-dependent stabilized jellium model can emulate 
quite well the second-order perturbative result of Rose and Dobson and also the non-perturbative 
calculation, with the single exception of the (110) case. 



4 Conclusions 

We have studied within perturbation theory surface energies of metallic slabs taking as 
zero-order the jellium model. Our results allow us to conclude that first-order surface energies 
depend strongly on the pseudopotential used. For aluminum slabs, the Heine-Abarenkov and the 
evanescent core potentials give similar results, while the Ashcroft potential differs from those two 
(it differs more for the least dense surfaces). It should therefore be used with some caution. On the 
other hand, the first-order perturbation turns out to be too large for the least dense surfaces, being 
imperative to correct it through second-order terms. We are implementing second-order perturba- 
tive calculations using the slab response function. We pointed out the importance of the Madelung 
energy, which should therefore be evaluated as far as possible without any approximations. 

We stressed the usefulness of the stabilized jellium model, which keeps the essential sim- 
plicity of jellium, while reproducing various first principles results for slabs and surfaces. As an 



extension of that model, we may take as zero's order the stabilized jellium energy instead of the 
jellium one (since the perturbation is strong for all aluminum faces, it should, after all, be a better 
starting point for perturbation theory) and to take as first-order perturbation, without averaging, 
the difference between the pseudopotential of stabilized jellium and an adequate pseudopotential. 
Work along these lines is in progress. 
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Figure captions 



Fig. 1: 

Three pseudo-potentials for aluminum used in this work as function of radial distance. 



Fig. 2: 

Jellium electronic density n (z) and perturbative potential Sv(z) contributing to first-order 
energy of slab with 7 layers (111 face). Three pictures refer, from top to bottom, to Ashcroft, 
Heine- Abarenkov and evanescent core pseudopotentials. The shaded area indicates the background 
jellium. 



Table 1: First-order surface energies for slabs with 7 and 19 layers of aluminum cut along, respec- 
tively, the (111), (100) and the (110) planes, aj, om, crps , and <tm,+j are the terms of the surface 
energy, a, following (2.15). The Perdew-Wang J2E] correlation functional was used in the jellium 
term. The other terms were calculated using the approximations expressed by (2.16), (2.17), and 
(2.18). All values are in erg/cm 2 . 



Face 


7 layers 
19 layers 


<?M 


Pseudopotential 


a {1) 

Vps 

7 layers 
19 layers 


%,+j 


a 
7 layers 
19 layers 


Al(lll) 


-610.4 
-602.7 


408.6 
408.6 


Ashcroft 


1050.5 
1053.5 






851.7 
856.4 


-610.4 
-602.7 


408.6 
408.6 


Heine- Abarenkov 


680.9 
677.8 






479.1 
483.7 


-610.4 
-602.7 


408.6 
408.6 


Evanescent core 


544.9 
543.5 


-92.7 
-92.7 


250.4 
256.7 


Al(100) 


-613.8 
-603.4 


1802.8 
1802.8 


Ashcroft 


405.3 
410.8 






1594.3 
1610.2 


-613.8 
-603.4 


1802.8 
1802.8 


Heine- Abarenkov 


-108.3 
-107.7 






1080.7 
1091.6 


-613.8 
-603.4 


1802.8 
1802.8 


Evanescent core 


-168.1 
-168.2 


-182.7 
-182.7 


838.1 
848.5 


Al(110) 


-585.6 
-608.3 


5540.3 
5540.3 


Ashcroft 


-1588.6 
-1581.7 






3366.1 
3350.3 


-585.6 
-608.3 


5540.3 
5540.3 


Heine- Abarenkov 


-2315.2 
-2306.2 


-15.4 
-15.4 


2624.1 
2610.4 


-585.6 
-608.3 


5540.3 

5540.3 


Evanescent core 


-1986.8 
-1978.1 


-588.7 
-588.7 


2379.2 
2365.1 



Table 2: Flat stabilized jellium surface energies (no corrugation included in the calculations) 
for slabs with thicknesses corresponding to 7 and 19 layers. <tsj,j arid <Jsj )PS are the terms of 
the stabilized jellium surface energy, asj, following (J2.20J) . The Perdew-Wang J2H| correlation 
functional was used in the jellium term. All values are in erg/cm 2 . 



thicknesss 


&SJ,J 


(?SJ,ps 


&SJ 


7 layers of Al(lll) 
19 layers of Al(l 11) 


-447.9 
-450.6 


1369.8 
1376.0 


921.9 
925.4 


7 layers of Al(100) 
19 layers of Al(100) 


-445.9 
-450.2 


1363.2 
1376.2 


917.3 
926.0 


7 layers of Al(llO) 
19 layers of Al(l 10) 


-460.1 
-451.5 


1405.1 

1370.8 


945.0 
919.3 



Table 3: Surface energies obtained by first-order perturbative theory (Lang-Kohn approach pQ 
with the Ashcroft pseudopotential), in the second column, and by other methods. The third 
column refers to flat stabilized jellium [Zj, a S j , the fourth to face-dependent stabilized jellium [TT] . 
<J s a f e , the fifth to the perturbative- variational method of Monnier and Perdew fSJ, (Jmp-, the sixth to 
the second-order perturbative theory of Rose and Dobson j^j, <Jrd, and the seventh to Schochlin, 
Bohnen, and Ho [27] first-principles non-perturbative results, usbh- in these works the semi- 
infinite system is considered except in the last one, where a slab with 12 layers was taken. All 
calculations were done in the Local Density Approximation. All values are in erg/cm 2 . 



Surface 


a 


flat 


face 
°S.J 


&MP 


VRD 


&SBH 


Al(lll) 


842 


925 


938 


643 


1065 


939 


Al(100) 


1631 


925 


1087 


1460 


1160 


1081 


Al(110) 


3393 


925 


1679 


2870 


1700 


1090 



h 

43 



c 
o 



Ashcroft 

Heine- Abarenkov 
Evanescent core 




Radial distance (bohr) 



Figure 2 
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